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Abstract. Radio observations in the ammonia inversion lines of four comets, C/2001 A2 (LINEAR), 153P/Ikeya-Zhang, 
C/2001 Q4 (NEAT) and C/2002 T7 (LINEAR), were performed at the Effelsberg 100-m Radio Telescope during their respective 
close approaches to Earth. None of the four lowest energy metastable lines (J, K = J), J = 1-4, could be detected in these 
comets. We derive the following 3tr upper bounds on the NH3 production rate , and comparing to the corresponding water 
production rates, percentage NH 3 abundances relative to H 2 Q: g(NH 3 ) < 1.9xl0 26 s -1 (0.63%) for C/2001 A2 (LINEAR), 



<2(NH 3 ) < 2.7xl0 2 ' 



(0.13%) for C/2001 Q4 (NEAT), g(NH 3 ) < 2.3xl0 2 



(0.74%) for C/2002 T7 (LINEAR) and 



g(NH 3 ) < 6.3xl0 2b s~ J (0.63%) for Comet 153P/Ikeya-Zhang. At 0.74% or less, the ammonia-to-water ratios are factors of ~2 
below the value for C/1995 01 (Hale-Bopp) and lP/Halley, suggesting chemical diversity between comets. The 18-cm lines of 
OH were clearly detected in the two comets observed during the 2004 campaign, thereby validating the cometary ephemerides. 

Key words. Comets: individual: Comet C/2001 A2 (LINEAR); Comet 153P/Ikeya-Zhang; Comet C/2001 Q4 (NEAT); Comet 
C/2002 T7 (LINEAR); Radio line: solar system 



1. Introduction 

Attempts to detect the radio K-band lines of ammonia (NH3) 
were made during the recent apparitions of four comets with 
the 100-m Effelsberg Radio Telescope of the Max-Planck- 
Institut fiir Radioastronomie (MPIfR). Observations were per- 
formed on Comets C/2001 A2 (LINEAR), 153P/Tkeya-Zhang, 
C/2001 Q4 (NEAT) and C/2002 T7 (LINEAR). A preliminary 
account of the observations of comets C/2001 A2 and 153P 
was previously published by Bird et al. (2002). Ammonia had 
previously been detected at MPIfR in the comets C/1983 HI 
(IRAS-Araki-Alcock) (Altenhoff et al. 1983) and C/1995 Ol 
(Hale-Bopp) (Bird et al. 1997, 1999), and at the Green Bank 
Telescope in C/1996 B2 (Hyakutake) (Palmer et al. 1996). 
Column density predictions from model calculations were en- 
couraging for these recent comets due to their relatively favor- 
able viewing geometries. 

Ammonia is expected to be among the more abundant 
volatile parent molecular constituents of cometary nuclei. 
Chemical models predict an abundance relative to H2O ice of 
a few percent (Charnley & Rodgers 2002 1. Our observational 
knowledge about NH3, ironically, comes largely from optical 
spectra of NH2, its longer-lived dissociation product. A recent 
reanalysis of the NH2 spectra of many comets suggests a typi- 
cal NH 3 abundance of 0.5% (Kawakita & Watanabe 2002, us- 
ing the comet database of Fink & Hicks 1996 1. 

This is less than the estimates of 1.0-1.8% from the ra- 
dio detections of NH 3 in C/1995 Ol (Hale-Bopp) (Bird et 



al. [TW7I |2TjU21 Hirota et al. but consistent with the 

revised estimate of 0.6% in C/1996 B2 (Hyakutake) (Bird et 
al. 119991 . The NH2 data imply an ammonia abundance of 
0.75% for comet lP/Halley, i.e. only about half the in situ value 
of 1.5% NH3 measured during the flyby of the Giotto space- 
craft (Meier et al. H994l >. 



In the interstellar medium, for comparison, the best esti- 
mates we have for the NH3 abundance in ice are upper lim- 
its of 5% and 7% on the lines of sight towards W33A (Taban 
et al. 2003 1 and the massive protostars GL 989 and GL 2136 
(Dartois et al. 2002 1, respectively. 



Few radio spectra of NH3 have been obtained because of 
the short photodissociation lifetime in the interplanary medium 
near 1 AU (r ^ 6000 s). The diameter of the ammonia cloud 
around comets is typically only a few thousand kilometers 
and beam dilution can be significant. The radio line observa- 
tions at MPIfR were thus scheduled near the relatively close 
approaches of the comets to Earth, thereby insuring that the 
source region of emission would be larger than, or at least only 
slightly smaller than, the 40" MPIfR beam. Contrary to expec- 
tations, no NH3 lines could be detected with certainty during 
the MPIfR comet observation campaigns in 2001-2004. The 
nondetections in C/2001 A2, 153P, C/2001 Q4 and C/2002 T7 
suggest an underabundance of NH3 in these comets with re- 
spect to Comet C/1995 Ol (Hale-Bopp) and lP/Halley. 
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Table 1. Orbital Elements: MPIfR Comet Observation 
Campaign 2001-2004 



Comet 



C/2001 A2-B 
(LINEAR) 



153P 
(Ikeya-Zhang) 



C/2001 Q4 
(NEAT) 



C/2002 T7 
(LINEAR) 



Orb. element set " 
Asc. node (o) 
Arg. perih. a> (o) 
Inclination (o) 
Perih. dist. q (AU) 
Eccentricity e 
Perih. epoch 
Osc. epoch * 



2001-M48 
295.12564 
295.32853 
36.47538 
0.7790280 
0.9993455 
2001 May 24.52062 
2001 May 11.0 



2002-H05 
93.37048 
93.37048 
28.12159 
0.5070583 
0.9899505 



2004-J69 
210.27852 
1.20653 

99.64258 
0.9619575 
1.0007438 



2004-J05 
94.85882 
157.73671 
160.658324 
0.6145967 
1.0005157 



2002 Mar 18.97927 2004 May 15.96707 2004 Apr 23.06172 



2002 Mar 27.0 



2004 May 8.0 



2004 Jun 4.0 



" Minor Planet Electronic Circular (MPEC) designation 
b Epoch of osculating elements 



Table 2. Mean Observation Parameters: MPIfR Comet Observation Campaign 2001-2004 



Comet Af 1 elevation RA.(J2000)_DEC A r /SEC 2 ZSCE 3 d 4 A6 5 _d 

date (min) (o) (hrmn) (o /) (AU) (AU) (o) (o) (10 3 km) (10 3 km) AG 



C/2001 A2-B (LINEAR) in 2001 



07.22 Jul 


70 


43 


23 31 


05 37 


0.262 


C/2002 CI (Ikeya-Zhang) in 2002 








25.42 Apr 


53 


57 


20 50 


61 29 


0.409 


04.40 May 


83 


58 


18 19 


55 30 


0.410 


06.38 May 


17 


51 


17 55 


52 38 


0.416 


07.41 May 


45 


54 


17 43 


51 14 


0.421 


C/2001 Q4 (NEAT) in 2004 








08.73 May 


47 


29 


07 44 


-07 53 


0.327 


13.69 May 


90 


53 


08 23 


+ 13 54 


0.390 


17.65 May 


82 


65 


08 47 


+25 42 


0.473 


C/2002 T7 (LINEAR) in 2004 








08.40 May 


35 


35 


00 41 


-4 06 


0.534 


13.44 May 


63 


28 


01 43 


-10 06 


0.364 



1.130 


109 


58 


9.7 


7.9 


1.22 


0.974 


74 


83 


9.5 


12.3 


0.77 


1.124 


95 


63 


11.8 


12.3 


0.96 


1.157 


100 


59 


12.3 


12.5 


0.98 


1.174 


102 


57 


12.6 


12.7 


0.99 


0.970 


74 


87 


9.3 


9.9 


0.94 


0.963 


72 


86 


9.1 


11.7 


0.78 


0.962 


70 


82 


9.1 


14.3 


0.64 


0.699 


41 


109 


5.4 


16.1 


0.34 


0.755 


38 


125 


6.0 


11.0 


0.55 



2. Ammonia Line Observations 

The 18-26 GHz HEMT receiver system (primary focus) was 
used for this program. The lowest four inversion transitions of 
ammonia in its metastable states (J, K = J), J — 1-4, were 
observed simultaneously in split-mode, dividing the autocorre- 
lator into 4 bands with 2048 channels each. Raw spectra were 
recorded at a typical channel spacing of 0.25 km s _1 over the 
range of at least ±30 km s about the expected line frequency 
at zero velocity in the comet rest frame. Spectra were cali- 
brated using W3(OH) and system temperatures were generally 
between 35 and 60 K. Various observing modes were used: fre- 
quency switching (C/2001 A2), position switching every 30s 
with a throw of either 5' or 10' (C/2001 A2, 153P), and beam 
switching using the rotating horn with a beam throw of 2' and 
a Is cycle (C/2001 Q4, C/2002 T7). TableHJpresents the orbital 
elements used to track the comet position and velocity. 

A summary of the mean geometric parameters during the 
observations, including integration time At, geocentric distance 
A, and heliocentric distance r, is given in Table [2] This table 
also shows values for the estimated diameter of the NH3 cloud, 
d — 6v ■ t ■ r 2 , with r the heliocentric distance (normalized to 



1 AU), t the mean lifetime of an ammonia molecule in inter- 
planetary space, and 5v the velocity spread of the coma gas as- 
suming spherically symmetric outflow. The mean lifetime of an 
ammonia molecule in interplanetary space depends on the solar 
UV flux and is probably somewhat longer than average during 
the 2004 epoch near solar minimum. In the following we use 
values for the photodissociation lifetime of NH3 at 1 AU in the 
range 5400 s< r < 5800 s. Different values for the line width 
were taken for each individual comet: 1.4 km/s (C/2001 A2); 
1.8 km/s (153P); 1.7 km/s (C/2001 Q4); and 1.9 km/s (C/2002 
T7). These were selected on the basis of independent radio line 
observations of other molecules during the same epoch (Biver 
et al., in preparation). The cloud diameter d may be compared 
to the next column, the physical extent of the FWHM antenna 
beam at MPIfR (AO, with 9 = 40") at the comet. The ratio of 
these quantities, d/A9, is shown in the last column of Table 
12 As a rule, the detection probablility is more favorable for 
large values of this ratio. The best conditions for this criterion 
held for comet C/2001 A2 and were still fairly good for comet 
153P. The beam was distinctly larger than the ammonia cloud 
for most of the observations in 2004. 



2001/A2 



Hatchell et al.: Searches for radio lines of NH3 
153P 2001/Q4 



2002/T7 



~i — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 



1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 



~i — 1 — 1 — 1 — 1 — 1 — r 



~i — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 



50 



(4,4) 






50 — 



(3,3) 

















1 1 


1 1 






1 1 


1 1 1 




50 



E-h 
50 



(2,2) 





(1.1) 




— 



-20 

v (km/s) 

Fig.l. NH 3 spectra of all comets (left to right: C/2001 A2, 153P, C/2001 Q4, C/2002 T7) in the lowest 4 metastable levels, 
centered at zero velocity in the comet rest frame. 



The Comet C/2001 A2 was observed to split into multiple 
nuclei (Sekanina et al., 2002). The position of the largest frag- 
ment (denoted 'A2-B'), targeted for these observations, was 
derived from the then current ephemeris (orbital elements in 
MPEC 2001-M48). The initial observations of Comet 153P 
on 25 Apr 2002, a long session in good weather, hinted at a 
marginal detection in the (1,1) and (3,3) lines (see Bird et al., 
2002). Follow-up observations performed on 4, 6 and 7 May 
under less favorable weather conditions (and higher system 
noise temperatures), however, could not confirm the tentative 
ammonia lines. Comet C/2001 Q4 could not be observed at 
MPIfR until early May 2004 due to its high southern declina- 
tion. The southward motion of Comet C/2002 T7 became prob- 
lematic as it approached the Earth, thereby severely restricting 
the last observation interval. In fact, no NH3 observations of 
Comet C/2002 T7 could be performed on 17 May 2004 due to 
lack of time following the OH observations (see next section). 



This and the unfavorable declination resulted in the shorter ac- 
cumulated integration time. 

The summed spectra for each ammonia line from all comet 
observations are displayed on the same scale (T MB in mK over 
±30 km/s) in Fig.^ The four spectra (from bottom to top) cor- 
respond to the (1,1) to (4,4) metastable states, and are centered 
at zero velocity in the comet rest frame. 

Upper limits to the line amplitudes are estimated from the 
spectra by fitting a Gaussian line model. This is preferable to 
taking the 3<x upper limit from the RMS noise level calculated 
over the whole spectrum because the actual data values in the 
central channels are taken into account. Therefore if there is 
any evidence for a weak line, that will increase our upper lim- 
its on the line amplitudes accordingly, in order to definitely ex- 
clude the weak line. This was the case for Comet 153P/ Ikeya- 
Zhang. 

The model spectrum consists of a Gaussian line of ampli- 
tude A plus a baseline level B. The baseline level is considered 
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Table 3. Ammonia line observations of comets at MPIfR and 
upper limits on the NH3 column density. 



IjT, 
[mK-km s 1 ] 



line 

(JJ) 



frequency 



T 

1 rms 

[mK] 



7*99.7% 

[mK] 



MB dV\ max 1 



(N(J,J)) max < 
[10 12 cm" 2 ] 



n(J,Jf (N(^)) max 
[10 12 cm" 2 ] 



Comet C/2001 A2-B (LINEAR): 7 Jul 2001 (At = 70 min, Sv ■■ 

(1,1) 
(2,2) 
(3,3) 
(4,4) 
Lowest 1 



16 
16 
16 
15 



Comet 

(1,1) 
(2,2) 
(3,3) 
(4,4) 
Lowest 

Comet 

(1,1) 
(2,2) 
(3,3) 
(4,4) 
Lowest 

Comet 

(1,1) 
(2,2) 
(3,3) 
(4,4) 
Lowest 



23.69450 
23.72263 
23.87013 
24.13942 
upper bound d 
153P/Ikeya-Zhang: 25 Apr 2002 ' 



36 
30 
34 
30 



54 
45 
51 
45 



23.69450 12 40 

23.72263 15 26 

23.87013 16 40 

24.13942 14 14 
upper bound d 

C/2001 Q4 (NEAT) 
23.69450 18 
23.72263 19 
23.87013 19 
24.13942 23 
upper bound d 

C/2002 T7 (LINEAR): 8,13 May 2004 (At 



(At = 53 min, Sv 
77 
50 
77 
27 



1.4 km s" 1 ) 
0.73 
0.46 
0.46 
0.38 

1.8 km s" 1 ) 
1.04 
0.51 
0.69 
0.23 



13,17 May 2004 (Af: 
20 36 
28 51 
44 80 
20 36 



219 min, Sv 



■ 98 min, Sv ■■ 



23.69450 
23.72263 
23.87013 
24.13942 
upper bound d 



29 
29 
29 
47 



30 
82 
52 
80 



61 
166 
105 
162 



1.7 km s" 1 ) 
0.49 
0.52 
0.72 
0.31 

1.9 km s" 1 ) 

0.82 

1.69 

0.95 

1.37 



0.26 
0.18 
0.20 
0.05 



0.26 
0.18 
0.20 
0.05 



0.26 
0.18 
0.20 
0.05 



0.26 
0.18 
0.20 
0.05 



< 2.8 

< 2.5 

< 2.3 
<7.6 

< 2.3 

< 4.0 

< 2.8 

< 3.5 

< 4.6 

< 2.8 

< 1.9 

< 2.9 

< 3.6 

< 6.2 

< 1.9 

< 3.2 

< 9.4 

< 4.8 

< 27.4 

< 3.2 



1.064 • 799.7% ' t> v with Sv values as given in the table. 

* from Eq. (TJ 

c defined in Eq. J5J 

d Lowest upper limit on (Af(NILi)) from the four lines observed 

e Calculations based on the best single epoch of observations 



again here as it is critical to the final line amplitudes, although 
linear baselines have already been subtracted from the spectra. 
The linewidth and line centre are fixed: the line centre always 
at km s _1 , and the linewidth taken from the observed radio 
linewidths of other molecules (Biver et al., in preparation; see 
Table [3] A grid in parameter space is calculated for the line 
amplitude A and the baseline level B. We calculate the proba- 
bility P(A, B) that the line has an amplitude A with baseline B 
using a likelihood function exp(~x 2 /2) with the constraint that 
A must be positive. As we are not interested in any residual 
baseline B, we integrate over B for each value of A to calcu- 
late the probability of each amplitude independent of baseline, 
P(A) = jP(A,B)dB a YjBj P(A t , Bj). From P(A) we identify 
the 99.7% upper limit on A (equivalent to 3cr for a normal dis- 
tribution). The corresponding upper limit on the line integrated 
intensity is calculated from 1.064 • A ■ Sv. A summary of the 
NH3 observations and derived results for all comets is given in 
Tabled 

3. Hydroxyl Line Observations 

Following the nondetection of NH3 lines during the observing 
sessions of Comets C/2001 Q4 and C/2002 T7 on 8 and 13 



May 2004, it was decided to verify the telescope comet track- 
ing configuration and, at the same time, the accuracy of the 
ephemerides. An imprompu change of receiver was requested 
and the initial hours of the final allocated observation interval 
on 17 May were devoted to observing the two strongest hyper- 
fine transitions of the A-doublet ground-state of OH. The OH 
spectra recorded for each comet are shown in Fig.|2] The upper 
(lower) panels show the two strongest OH lines at 1667 (1665) 
MHz, centered at zero velocity in the rest frame of the comet 
ephemeris. The total integration times, for comparison with the 
NH3 observations, were only 12 and 16 minutes on the Comets 
C/2001 Q4 and C/2002 T7, respectively. 

A tabular summary of these observations, which were the 
first cometary OH-observations performed at MPIfR since 
Comet Halley (Bird et al. 1986), is given in Table H The in- 
tensity of the A-doublet lines is governed to a large extent by 
the inversion parameter, a measure of the imbalance in the up- 
per and lower levels of the OH ground state (Schleicher and 
A'Hearn, 1988). If the two levels are nearly equal, the inversion 
parameter is zero and the OH-maser is deactivated. Although 
this was very nearly the case for both comets on 17 May 2004, 
clear detections were obtained for all combinations except for 
comet C/2001 Q4 at 1665 MHz. The 1667 MHz line is about 
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Table 4. OH Line Observations of Comets at MPIfR, 17 May 
2004 



transition 
F„ -> F, 



frequency 
[MHz] 



* peak 

[mK] 



Vo 

[km s" 1 ] 



FWHM = 6v 
[km s 1 ] 



Line Area " 
[mJy-km s _1 ] 



Comet C/2001 Q4 (NEAT): 17 May 2004 (Af = 12 min) 
1->1 1665.4018 72* 

2^ 2 1667.3590 119 -0.16±0.26 2.06+0.55 
Comet C/2002 T7 (LINEAR): 17 May 2004 (Af = 16 min) 
1 -> 1 1665.4018 99 0.54±0.23 2.69+0.36 
2^ 2 1667.3590 175 0.44+0.22 2.18+0.47 



130+31 

142+23 
203+39 



" 7V 



dv (MPIfR gain: 2.0 K/Jy at 1665/1667 MHz) 



2002/Q4 2002/T7 



1667 




-20 20 

v (km/s) 

Fig. 2. OH radio line spectra recorded on 17 May 2004. 
Left: C/2001 Q4 (NEAT), 12 minutes integration time; Right: 
C/2002 T7 (LINEAR), 16 minutes integration time. The main 
beam brightness temperature T MB is shown for the two OH lines 
at 1667 MHz (upper) and 1665 MHz (lower). 



twice as strong as the 1665 MHz line in both comets, consistent 
with the 9:5 ratio expected under LTE conditions. 

It was verified that the OH 1667 MHz line strengths were 
consistent with simultaneous measurements of the radio OH 
lines taken at the Nancay Radio Telescope (J. Crovisier, pri- 
vate communication, see Table [5}. Note that the Nangay line 
strengths given in the table represent an average over the 1665 
and 1667 MHz lines, both polarizations (LCP+RCP), and then 
scaled to the 1667 MHz line, assuming they conform to the 
LTE ratio. The line strengths at MPIfR were found to be some- 
what higher than at Nangay, but are within the l<x measure- 
ment errors. An exact agreement would not be expected as al- 
though the Nangay and MPIfR beams at 18cm are similar in 
area, they are very different shapes and therefore couple to the 
comet OH emission differently (the Nangay beam is elongated 
with FWHP 3.5' x 19' , whereas the MPIfR beam is approxi- 
mately round with a FWHP 7.8'). 



Table 5. OH Comet Observations, Nancay RT and MPIfR, May 
2004 



May Date Line Area ° 


Line Area * 


Line Area c 


LTE 


LTE 


1667 only 


[mJy-km s ] 


[mJy-km s _1 ] 


[mJy-km s -1 ] 


Comet C/2001 Q4 (NEAT) 






16 104+12 


no obs. 


no obs. 


17 101+12 


65-207 


130+31 


18 75+12 


no obs. 


no obs. 


Comet C/2002 T7 (LINEAR) 






16 176+12 


no obs. 


no obs. 


17 no obs. 


229+28 


203+39 


18 185+16 


no obs. 


no obs. 



" Nancay RT, mean of 1665 MHz and 1667 MHz lines assuming 
LTE ratio of 5:9 

h MPIfR, mean of 1665 MHz and 1667 MHz lines assuming LTE 
ratio of 5:9. Limits for C/2001 Q4 are assuming the 1665 MHz line 
amplitude is (a) (b) 3cr upper limit. 

' MPIfR, 1667 MHz detection only 



4. Upper limits on NH 3 column densities and 
production rates 

Table 3 (column 6) shows the upper bounds on the column 
density {N{ J, J)) ma x (in cirT 2 ) of molecules in the given state. 
These are calculated from the line integrated intensity in Kkm 
s _1 using the expression (e.g., Rohlfs & Wilson 1996, p. 362) 



(N(J,J)) = 6.8x10' 



,J+ 1 
J 



T MB (v)dv 



(1) 



Both the upper and lower levels of the specific metastable NH3 
state are included in Eq. Q. 

Knowing the mean column density of the NH3 molecules 
in a given metastable state, the total column density of all NH3 
molecules is obtained from the relation 



n(J,J) = (N(J,J))/(N(NH i )) 



(2) 



where n( J, J), the relative population for each observed state, 
depends on the mean kinetic temperature and density of the 
cometary NH3 gas in the antenna beam. The population will 
favor the metastable (K=J) states, particularly if the collision 
time is much longer than the decay time scales ~10 s for the 
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Table 6. Comet Production Model Parameters 



Comet 


mm vvvv 
j j j j 


Qo 


r 


Sv 




T 


Reference 






[10 29 s" 1 ] 


[AU] 


[km s" 1 ] 


[s] 






C/1995 01 (Hale-Bopp) 


04 1997 


100 


0.914 


1.8 


3.0 


5600 


Colometal. (1999) 


C/2001 A2-B (LINEAR) 


07 2001 


0.28 


1.130 


1.49 


3.0 


5400 


E Bensch (priv. comm.) 


153P/Ikeya-Zhang 


04 2002 


0.92 


1.000 


1.77 


3.21 


5600 


Dello Russo et al. (2004) 


C/2001 Q4 (NEAT) 


05 2004 


1.8 


1.017 


1.67 


3.0 


5800 


N. Biver (priv. comm.) 


C/2002 T7 (LINEAR) 


05 2004 


4.5 


0.645 


1.59 


3.0 


5800 


N. Biver (priv. comm.) 



nonmetastable states. We have performed statistical equilib- 
rium calculations using collision rate coefficients of ammo- 
nia in an H2 background gas (Danby et al. 1988 Schoier et 
al. 12005ft . but multiplied by a factor 4.3 to account roughly 
for the increase in cross section for collisions with H2O. The 
larger mass of H2O than H2, which decreases the collision ve- 
locity, was not taken into account. The calculations included 
the lowest 17 states for ortho NH 3 up to (J.K)= (6,0), 599 K 
above ground, and the lowest 24 states for para NH3 up to 
(5,1), 423 K above ground. The population of higher energy 
states was found to be negligible for densities up to 10 9 cm" 3 
and kinetic temperatures up to 200 K. The ortho-to-para ratio 
was assumed to be unity, but may, in fact, be slightly higher 
(Kawakita et al., 2001). The coma gas density for the observed 
comets, based on their later determined production rates, lie 
in the range from 10 5 to 10 6 cirT 3 at a distance R = 5000 km 
from the nucleus. This can be taken as a typical mean den- 
sity of the background gas in the telescope beam. The model 
calculations show that the NH3 partition function in the coma 
is fairly insensitive to the actual density and only moderately 
dependent on the kinetic temperature over the range of plausi- 
ble values from 50 K < < 100 K. The calculated values of 
n(J, J) in statistical equilibrium with 7j = 100 K are shown in 
Table |3 for each transition. The remaining column of Table [3] 
presents the total NH3 column densities, (N(NHi,)) max , calcu- 
lated from Eq. for the upper bounds of each transition and 
comet given in Table [3] To calculate the production rates, we 
use the strongest upper limits on /V(NH3) from among the four 
individual NH3 lines. 

The upper bounds on the column densities for all observed 
comets are plotted in Fig. [3] and compared with pre-perihelion 
model predictions of (N(NHt,)) over an interval +90 days 
about perihelion (curves). Fig.|3]also shows the measurements 
and prediction curve for the firm detection of NH3 in Comet 
C/1995 Ol (Hale-Bopp) (solid black line and filled circles). 

The model calculations assume that the ammonia produc- 
tion rate follows the water production rate with an abundance 
ratio of exactly 1 .0%, a mean lifetime of t « 5600 s (with small 
corrections for the phase of the solar cycle), and an antenna 
beamwidth of = 41.5". 

The H2O production rate and gas outflow velocity for each 
comet were assumed to vary with heliocentric distance r (in 
AU) according to (e.g., A'Hearn et al. 1995): 



GW = Qo 



6v{r) = 



(3) 



the linewidths at each observation epoch (6v) given in Sec. [2] 
The model parameters and the reference sources used to calcu- 
late the prediction curves of Fig. [3] are listed in Table [6] 

The unique stature of C/1995 Ol (Hale-Bopp), which had 
a production rate much greater than the other comets, is clearly 
evident in Fig. [5] and Table [6] The solid circle points with rms 
error bars for C/1995 Ol (Hale-Bopp) in Fig. [3] are derived 
from the measured (3,3) line strengths (Bird et al. 1997). The 
water production rate at perihelion, 1.0 x 10 31 s _1 , was taken 
fromColom et al. (1999). 

The procedure originally developed by Snyder (1982) is 
used here to estimate the production rate Q(S ,) of species i 
from the observed beam-averaged column density (N(S ,)) 



(N(Sd) = 



4-Q(Sj) ( d/A0 for d < AO 
ndv Ae{T(d/A0) for d > AO 



(4) 



with d the comet coma diameter, A the geocentric distance, 6v 
the linewidth (taken as 2m, where u is the gas outflow velocity), 
and the half-power beam width of the antenna. The function 
T , given by 



!M\ d 

T = arccos — H 

1 d / A0 




(5) 



with a (usually a » 3) and the values of r„, and Q„ taken from 
the literature. The linewidths at 1 AU (6v n ) are converted from 



varies from T - 1 for d = AO to T = n/2 for d » AO. 

Using Eq. @, we can derive an upper bound on a comet's 
NH 3 production rate from the upper bounds on the column den- 
sities given in the bottom lines for each comet of Table [3] The 
resulting NH3 production rates for these comets are given in 
Tableland range between 1.8 and 23 x 10 26 s" 1 . 

5. Ammonia-to-water ratio 

The strongly variable water production rate for Comet 
C/2001 A2-B was reported as g(H 2 0) = 3.8 x 10 28 s" 1 
on 1.7-2.0 July 2001 (Biver et al. 2001), consistent also 
with g(H 2 0) ~ 4.0 x 10 28 s _I observed on 9-10 July 2001 
(Dello Russo et al. 2005). The Odin satellite measured 5- 
8 x 10 28 molec. s" 1 between 20 June 2001 and 7 July 2001 
(Lecacheux et al. 2003). A recent compilation of contin- 
uous water production rates derived from SWAS observa- 
tions (Submillimeter Wave Astronomy Satellite) of Comet 
C/2001 A2-B (F. Bensch, private communication) suggests a 
value of 2(H 2 0) = 2.8 x 10 28 s" 1 as appropriate for the MPIfR 
observations on 7 July 2001. These water production rates im- 
ply that the maximum NH3 fraction for the Comet C/2001 A2- 
B observations becomes 0.63%. 
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Fig. 3. Predicted and observed MPIfR beam-averaged NH3 column densities for five comets. The thick part of the curves denote 
intervals where the declination enables observation at MPIfR (5 > -25°). Maxima occur near closest approach to Earth. Points 
with standard deviations represent the line detections in C/1995 Ol (Hale-Bopp) (Bird et al. 1997); points without error bars are 
nondetection upper bounds for the Comets C/2001 A2-B, 153P, C/2001 Q4 and C/2002 T7. In all nondetection cases, the upper 
bounds are significantly less than the predictions based on an ammonia abundance of 1%. The shortfall is indicated in each case 
by the dotted line. 



An estimate of the water production rate in Comet 
153P/Ikeya-Zhang of g(H 2 0) = 1.7 x 10 29 s 1 could be de- 
rived from Odin satellite observations of the 557 GHz line of 
water on 26.8 April 2002 (Crovisier et al., private commu- 
nication; Lecacheux et al. 2003). HST ultraviolet OH obser- 
vations on 20-22 April 2002 (Weaver et al., private commu- 
nication) yielded £>(H 2 0) = 2.3 x 10 29 s _1 . Newer system- 
atic estimates of the H2O production for Comet 153/P Ikeya- 
Zhang (Dello Russo et al. 2004) are smaller by more than a 
factor of two with respect to the original estimates used in the 
NH3 abundance estimates of Bird et al. (2002). The best fit 
model prediction of Dello Russo et al. (2004), the parameters 
of which are given in Table [6] yields a water production rate of 
g(H 2 0) = 1.0 X 10 28 s 1 for the mean epoch of the MPIfR ob- 
servations. We thus calculate an upper bound on the ammonia- 
to-water ratio in Comet 153/P Ikeya-Zhang of 0.63%. 

Water production rates for the Comets C/2001 Q4 and 
C/2002 T7 have been reported (N. Biver, private communica- 
tion) as 0(H 2 O) = 1.8 + 0.2 x 10 29 s" 1 on 27.0 April 2004, 
and 2(H 2 0) = 4.5 ± 1.0 x 10 29 s" 1 on 2.0 May 2004, respec- 
tively. These lead to upper bounds on the ammonia fractional 
abundance relative to water of, respectively, 0.13% and 0.74%. 



The upper limits for the relative NH3 abundance in each 
comet, based on these H 2 production rates and the NH3 pro- 
duction rates derived in Sec.|4l are presented in Table For 
comparison, a summary of the ammonia abundance estimates 
from all previous comet observations has been compiled by 
Bockelee-Morvan et al. (2005). 

All of the upper limits we present here are significantly 
lower than the estimate of 1.1% for C/1995 Ol (Hale-Bopp) 
(Bird et al. 1999) or 1.5% for Comet lP/Halley (Meier et 
al. 1994). The relative abundance implied for C/2001 Q4 was 
found to be considerably lower than the values of 0.64-0.74% 
derived for the other three comets. This may be a real effect, 
but is subject to reconfirmation when finally agreed values for 
the water production rate become available. 

6. Conclusions 

Our observations suggest that there is a diversity of the NH3 
abundance among comets. This was already suggested by the 
large (more than a factor of ten) range of the NH/OH ratio ob- 
served from narrowband photometry in the visible for a large 
sample of comets by A'Hearn et al. (1995). The nondetections 
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Table 7. NH3 production rates and abundances relative to H2O 



Comet 


epoch 


g(NH 3 ) 


e(H 2 o) r 


[NH 3 ]/[H 2 0] 






[10 26 s" 1 ] 


[10 29 s- 1 ] 




C/2001 A2-B (LINEAR) 


7.22 Jul 


1.8 


0.28 


< 0.63% 


153P/Ikeya-Zhang 


25.42 Apr 


6.4 


1.0 


< 0.63% 


C/2001 Q4 (NEAT) 


14.11 May 


2.7 


2.1 


< 0.13% 


C/2002 T7 (LINEAR) 


11.64 May 


23 


3.1 


< 0.74% 



of ammonia in C/2001 A2-B (LINEAR), C/2001 Q4 (NEAT) 
and C/2002 T7 (LINEAR), as well as the very marginal detec- 
tion in P153 Ikeya-Zhang, lead to the preliminary indication 
that the NH3 fraction may be a factor of ~2 lower than that de- 
rived from the detections in Comets lP/Halley and C/1995 Ol 
(Hale-Bopp). 

With the possible exception of C/2001 Q4 (NEAT), it is in- 
teresting that the upper limits all seem to be consistent with 
the 0.5% abundance determined from the NH2 analysis of 
Kawakita & Watanabe (2002). 

The accumulated statistical sample is clearly insufficient 
for drawing a final conclusion on the spread in relative abun- 
dances of NH3 in cometary comae. Nevertheless, the current 
comet count does imply that the amount of ammonia in the 
archetypical comets lP/Halley and C/1995 Ol (Hale-Bopp) 
may be the exception rather than the rule. Only more obser- 
vations will provide final resolution of the issue. 
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